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Abstract 



This report examines design considerations of a flow system for 
transient electrical laser experiments. After a brief review of corona 
discharge gas flow interaction phenomena, the choice of the size and 
operating conditions of the discharge cell are discussed. Once chose 
parameters are established, the flow system, which consists essentially 
of a Ludwieg tube, is designed. Attention is paid to flow starting 
times, to the available time of steady state flow, and to the effects 
of the boundary layer on the core flow. Suggested experiments and in- 
strumentation are given, as well, in this report. 
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I . INTRODUCTION 



The effect of gasdynamic means for stabilizing corona discharges 
in atmospheric air has been shown ^'^#3) icantly important. 

By choosing the right turbulence spectrum and flow speed, the electri- 
cal power fed into a diffuse corona discharge has been improved by as 
much as 500 times above the maximum no-flow power. The influence of 
turbulence, which is more substantial at the higher pressures (~1 atm) , 
is to promote better mixing and homogeneity of the gas. Thus, by gas- 
dynamic means, molecular gases are able to accept higher diffuse dis- 
charge power, making such discharge attractive for electric gas laser 
operation because of the promise of attractive output power capabili- 
ties. 

The temporal nonuniformities of the gas density created by the 
turbulence may cause phase distortions in the laser output beam, and 
thus deteriorate the laser performance, (due to the effect of the phase 
distortions on the far field intensity distribution) . However, turbu- 
lent corona discharge studies show that the strongly decaying turbulence 

(2 4) 

mainly affects the discharge at the region close to anode tips ' 

where the electric field is strong and the plasma is formed; this is 

the corona region of length L , and has a strong nonuniform electric 

a 

field E (Fig. 1) . The external zone is a quasi-neutral zone which has 
a 

a length of 1^, 1^ » L , and an electric field E^ which is considera- 
bly lower and more uniform than E . Thus, by exciting the turbulence 

cl 

at the anode tips region and letting the turbulence decay and become 
very weak at the external zone, where laser action occurs, the density 
non uni formity problem can be reduced to minimum. 
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The separation of the discharge region into two main regions with space 

scales L and (where << L^) and electric fields and E^ (where 

E >:> E^) appears to be very useful for efficient laser operation. 

Those regions may be controlled, almost independently, by changing 

electrode geometry (anode tips radii, electrode separation, electrode 

material, etc,) and flow characteristics. It is possible to choose the 

average electric field in the external region so that the E/N value 

will be in such a range that most of the electron energy will go to 

excite the upper laser level; as for example, in a CO ^ laser containing 

a mixture of He/^/cc^ a volumetric ratio of 3/2/1, more than 80% of 

the electron energy goes into the (001) and N 2 (v=l) vibrational 

levels for the value of E/N - 2 x 10 ^ v cm 2 in this respect, 

an analogy can be made between a c.w. electron -be am- sustained electric 

discharge and a c.w. gas dynamic stabilized corona discharge; both cases 

(7) 

operate in the Ionizer/Sustainer mode , which means that there is a 
separation between the electric field that produces the electrons and 
the sustained discharge, where the electrons acquire their optimal 
temperature Te (Te is determined by E/N) for laser action. 

In this report we propose an experimental study of the effects of 
gasdynamic means (mainly turbulence and flow speed) and different geo- 
metrical parameters of the electrodes on corona discharges for laser 
applications. 

The experimental system proposed is a versatile inexpensive set-up, 
in which the following can be studied. 

Fluid dynamics measurements ; 

Turbulence spectrum Flow rates 

Temperatures 



Velocity profiles Pressures 
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Electrical measurements; 



V-I characteristics 

probes for E, Ne , and measurements 
Optical measurements : 

Photography 
Schlieren photography 
In terf e rograms 
Light emission spectroscopy 
Laser performance measurements : 

Small signal gain 
Output power measurements 
Parameters to be studied : 

Turbulence intensity and spectrum 
Flow rates (speed and pressure) 

Gas mixtures 

Discharge electrical energy 
Distance between electrodes 
Electrodes geometry 

Preionizing devices , like u.v. radiation 

To generate the flow through the corona discharge , the application 

( 8 ) 

of the Ludwieg tube principle has been considered to be an attractive 

proposition because of its simplicity , low cost and also high flow 

quality. The useful flow time of this device depends mainly on the 

charge tube length and on the speed of sound of the gas in the tube (see 

mse c 

Section III). For air this time is approximately 5 — — . For a 

** 1 meter length 

charge tube of the order of 4m in length , the net useful flow time will 
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be 15-20 msec, which corresponds to -30 cavity flow times (for 5 cm 

cavity length and gas speed of 100 m/sec) . This procedure effectively 

simulates cw operation without requiring the heat sink electrodes / 

large compressors , heat exchanger etc./ that would be necessary for 

true c.w. operation. All other relevant characteristic times are much 

( 2 ) 

shorter than the flow time 

Section II deals with the basic design and sizing of the test cell. 
Once the cell dimensions/ the pressure of the gas and its flow speed are 
determined/ the Ludwieg tube flow system can be designed. The detailed 
considerations of the design are given in Section III. The overall ex- 
perimental set-up is shown schematically in Fig. 2. 
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II. TEST CELL DESIGN 



The test section design is basically the same as that of the cell 

(2 3 ) 

used previously in the gas dynamic -discharge interaction studies 9 
Some changes have to be made due to the use of this cell for laser ex- 
periments and as being a part of a different flow system. 

1. Cell Operating Conditions 

Gas mixtures Different gas mixtures are going to be used in 
this system: N^ / He, CO^ , ••• etc * As f ar as a system 

design is concerned, the difference among those mixtures is their 
specific heat ratio y and molecular weight* 

will consider two values for y: 
y = 1.4 (diatomic molecular gas) 
y = 1.67 (monatomic gas) 

Gas pressure The cell will be operated at pressures around 500 
torr which is the pressure at which maximum electrical power can be 
dumped into the corona discharge in no flowing air wLfch 4cm elec- 
trode gap^ . For different gases, different electrode gaps and for 
flowing gases the optimum operating pressure varies, thus the sys- 
tem will be designed to operate in the pressure range of 200 <_ p < 

800 Torr (lower and higher pressures can be achieved as well) . 

Flow speed The gas dynamics -discharge interaction becomes effec- 
tive at flow speeds of ^ 100 m/sec^^ . The faster the flow is, the more 
discharge performance improves. The cell will be designed to have a 
flow of about Mach number M = 0.5. 
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2. Cell dimensions 



Cross section dimensions: 

Cell height: d = 2.22” , 5.5 cm . 

Cell width: b = 6.5" ,16 cm. 

Dimension along the flow direction: 

Cell length: c = 4" , 10 cm , 

Dimensions d and c are the same as those in Ref. 3; however, 

the cell width, b, is bigger since along this direction the small 

signal gain (s.s.g.) is measured and laser cavity will be aligned. 

(6 9 ) 

The s.s.g. reported in the literature v ' for electrical lasers, 
range from 1% to 4% cm 16 cm cavity length should be adequate for 
gain and power measurements. 

3. Nozzle throat 

The cross section area of the discharge cell is S = b x d = 

2 

88 cm . In order to obtaxn M = 0.5 in the test section, the cross 
section area A*, at the throat for the different values of y, are 
given in Table I: 





y = 1.4 


Y = 1.67 


A/A* 


1.34 


1.31 


2, 

A* (cm ) 

.. 


66 

.. _ ... 


67 


h* (cm) 


4.1 


4.2 



Table I give also the values of throat height h*, for a one 
dimensional nozzle. 

Turbulence generating plates have to be inserted ahead of 
the electrode tips. Since the effective cross section of those 
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plates is smaller than that of the cell, their effective cross 
section should exceed the values of A* given in Table I; other- 
wise, the flow will choke at the plate, prior to reaching the test 
section. 

The turbulent plate used in previous experiments had an 
area blockage of the order of 50%, which is much above the 34% 
allowed for this design* A new design and location of those 
plates have to be considered. One way of improving the plate 
blockage is to mount the turbulence generators closer to the 
electrode tips. This will enable a weaker turbulence generation 
, (thus, a smaller blockage) since the decay of the turbulence at 
the electrode tips will be insignificantly small. One way of 
doing it is shown in Fig. 3. 

4* Electrical power heat addition 

The effect of adding heat to a compressible flow is to send 

compression waves into the flow which change flow parameters; heat 

addition to subsonic flow will accelerate the flow and for a cri- 

( 10 ) 

tical value of heat, q , the flow will be choked. 

tr 

o = q ° r = 1 (l + YM 2 ) 2 _ 

p stag 2(y + 1)M 1 + "" M 
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Typical mass flow rates through the cell: 

P 



m = pAv = 



(R/M) T 



Av 



±2. 4 

for N s A = 2 x 88 x 2 x L ? - = l Kg/sec 

( 3,3 2 q 10 - ■)« 300 



i£L 4 

for H : m = 2 x 88 .1 ° . ■ = 0 .14 Kg/sec 

6 ( e - J. 4 ,iQ .. ) K300 



Typical critical heat addition values are given in Table II: 





Q 

*cr 


q kw/. . 
^cr kg/sec 


q , kw } For typical 
m values 


for N 2 


0.45 


140 


140 


for B 

e 


0.4' 


620 


87 



Table II 



The power supply available has - 1 kw power output which is 
very low power compared to q . 
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III. LUCMIEG TUBE 



1. Introduction 

The Ludwieg tube flow system is shown schematically in Figure 2. 

The Ludwieg tube consists of a long tube, charge tube, where the gas 
is initially compressed to the desired charge pressure, P^. A tran- 
sition section channels the flow from the circular charge tube to 
the rectangular test section which is a part of a converging-diverging 
supersonic nozzle. The gas expanding through the nozzle is discharged 
through a short diffuser into a dump tank. A diaphragm separates the 
driver and driven sections. The diaphragm can be located either down- 
stream or upstream of the nozzle. 

The flow is initiated by the rupture of the diaphragm. In the 
case of an upstream diaphragm, a backwards facing centered expansion 
fan propagates into the charge tube. This fan is reflected at the 
closed end of the tube and returns back to the nozzle. This is called 
the first flow cycle; a simplified wave diagram is shown in Fig. 4(a) . 
Between the time when the tail of the expansion wave leaves the nozzle 
entrance and when the head of the reflected fan reaches the nozzle 
again, the stagnation conditions ahead of the nozzle remain constant and 
a clean steady flow may be established in the nozzle. A second flow 
cycle (and others) follows the first one but the stagnation parameters 
are different. The duration of the steady flow, for a certain gas, is 
primarily a function of the charge tube length and we are interested to 
make it as long as possible; there is a limit on the length of the tube 
set by the growth of the boundary layer ^ . The boundary- layer 

growth sets also a limit on the diameter of the tube and on the flow Mach 
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number, , in the charge tube. Subsequent to the diaphragm rupture, 
a shock wave and a contact surface proceed downstream through the 
driven section. Additional waves follow the contact surface to adjust 
the nozzle test section pressure to the freestream pressure correspond- 
ing to the designed Mach number expected at the nozzle exit for a given 
nozzle . The steady supersonic flow in the nozzle is not established 
until all those waves are swept out the nozzle. The time required to 
get the flow started, that is, the time required’ to achieve a steady 



flow from rest, in the supersonic nozzle, is a very important performance 
criterion. It is important to minimize the start time in order to maxi- 
mize the duration of the steady state conditions for a fixed charge tube 

. _ (14,15) 

length 

For downstream diaphragm location, the wave diagram is different 

from the case where the diaphragm is located upstream, a simplified x - t 

diagram is shown in Fig. 4(b) . As soon as the local sonic velocity is 

reached at the nozzle throat, no further expansion waves can pass through 

into the tube and an expansion fan of finite width passes down the tube. 

For an upstream diaphragm location, the starting times at the 

nozzle exit are reported to be shorter than those obtained with the 

(15) 

downstream diaphragm location . However, the disadvantage of the up- 
stream location diaphragm is that when breaking the diaphragm, particles 
might tear off and disturb the flow in the test section. 

2. Charge tube diameter and length 

Economic considerations dictate that a minimum charge tube diameter 
be used for a given test section size. This results in a relatively 
small contraction (small ratio between nozzle throat area to charge tube 
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area) which requires high Mach number/ , at the charge tube* This 
will cause the boundary layer to have a dominant effect on nozzle flow* 
It is suggested to have the I.D* of the charge tube to be the same 
as the dimension b of the test cell, i.e. 16 cm* In this case A^/A* 
and are given in Table III. 





Y = 1.4 


Y = 1.67 


V A * 


3.05 


3.00 


M 3 


0*2 


0.19 



Table III 



If a tube with a larger diameter is available (like 20cm) , it will im- 
prove the performance of the system, since will be smaller* It is 
not recommended to use a plenum chamber upstream of the nozzle inlet 
(to slow and mix th.e core and boundary layer flows emerging from the 
pipe) since it will be shown that for 5 0.2 the disturbances intro- 
duced by the boundary layer growth are small while the starting time is 
going to be much longer by introducing the plenum chamber. 

5*/d is given as a function of instananeou^^spect ratio of the 

by Russel , based 



flow s/d, P^d and M , in the form of 5*/d 

,(16) 






(s/d) 



4/5 



on Becker's model for 6*/d , where <5* is the boundary layer displace- 

ment thickness, at a distance s from the expansion wave head, and d is 
the charge tube diameter* 



for M_ = 0.2 -r 
3 d 



5* <V> 



1/5 



4/5 



= 5 x 10 4 (atm - m) 



(s/d) 

At the end of the test time, s = 2L and the thickness of the boun- 



dary layer is 6 = 105*. For = 1 atm, L = 3 m and d = 0.16 m. 
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we get 5*/d = 1.3%, 6 ~ 0.02 m and 6 /d ~ 13% . 

The perturbations in the core flow due to the boundary layer, 
when the test period is over, are (for the same parameters used for 
6 * 

— calculations) : 
d 



_AP 3 (P 4 d) 



1/5 



= 4.3 x 10 " 4 ^ = 6.5 x 10‘ 



3 (L/d ) 4 / 5 



AP. 



S V (P 4 d > V5 _,, -3 * ^-4.8*10 

v 4/5 x ^ ^3 

3 (L/ d) 



-2 



AT Y - 1 AP 

— - -r = - Ta x 6.5 x 10 = -1.7 X 10 

T 3 y P 3 1.4 



which are relatively very small disturbances. 

The effects of the above described core perturbations on the down- 
stream flow were discussed by Russell . He evaluated the boundary 
layer along the walls of a set of specified nozzles. The results for 
different values of 6 */d, M 3 and M (Mach number in the nozzle) are 
presented in Fig. 7 of his paper; he presents the effective area - 
(A/A*) ef f — A/A* x 

as a function of &*/d for different values of M 3 and 

Mj where (A/A*)^£ Is the test section - to - throat area ratio with 
boundary-layer correction. 

A/A* is the geometric area ratio, with no corrections. 
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For the value of P.d = 0 o 5 atm-m, Russell found that with 
4 

-4 * 

10 o*/d <_ 1%, the effective area change parameter has a finite 

value determined by the normal steady boundary layer correction of 

the nozzle itself (in our case P d - 2 x 0.16 0.3. There is a 

o 

weak dependence on pressure level (P^d) . 

To summarize: 

for M = 0 .5 Mach number flow in the test section, d=16cm, L=3m 
and m 3 = 0.2 

6*/d = 1.3% , 6/d = 13% 




0 . 6 % 




(A/A*) eff - A/A* 
A/A* 



0 



4% 



(contribution from charge tube 
boundary layer) 



3 . Test time 

The test time, (the time taken for the head of the reflected 
expansion wave to reach back the throat) for small values of is 
given approximately by 2L/a^, where a^ is the sound speed of the 
undisturbed gas in the charge tube. 

The time available falls only slightly with increasing the Mach 

. M (17) 
number M^ 



t/t o = rr 2 sr {1 + ;L T JiM 3 )2l ’ , ' 1) - 



where t = L/a. . For M- = 0.2 and L = 3m, the calculated values of 
O 4 J 



t and t are given in Table IV. 
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Y = 1.4 


Y = 1.67 


a^ (m/sec) 


350 


1000 


t (msec) 
o 


8.55 


3 


t/fc o 


1.88 


1.9 


t (msec) 


16 


5.7 



Table IV 



4, Starting times 

The starting time of supersonic flow at a given location in the 
flow system is defined as the time elapsed between the arrival of the 
head of the expansion fan (downstream diaphragm) , or shock wave 
(upstream diaphragm) and the instance where the pressure becomes time in 
dependent, i.e., when steady flow is established. For an upstream 
diaphragm location, the starting times at the nozzle exit are reported 
to be shorter than those obtained with the downstream diaphragm loca- 



a. Downstream diaphragm location: 

According to ref. 15, (they used a Ludwieg 'Tube with a charge tube 
tube of 670 cm in length and I.D. of 13.4 cm) the starting times at the 
nozzle exit are independent of the initial downstream pressure p^, 
except in the range above performance limit. Starting times increase 
with diffuser length. The starting times upstream of the nozzle throat 
do not depend noticeably on diffuser length, and are of the order of 
0.5 msec, where the starting times in the diverging part of the nozzle 
are of the order of 6 msec. The experimental results are presented in 
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dimensionless starting time t g /t^ as a function of dimensionless axial 
distance x/d, where is defined as &/**; Z' = /R*h*, R* is the radius 
of curvature at the nozzle throat, h* the nozzle height and a* the 
speed of sound at the throat. At the subsonic nozzle section t /x^ - 3 anc ^ 
at the nozzle exit t /x^ » 40. In our case, the starting time in the 
subsonic section will be t “ 1 msec and in the supersonic section 

t ~ 8 msec, for For other gases, the starting time scales in- 

versely to a*. 

b. Upstream diaphragm location 

The starting times, at the exit of the nozzle, increase with 
initial downstream pressure, P^ , and increase also with diffuser length. 

The starting times are of the order of 2 msec at nozzle exit and 
should be at least of the order of 10 in order to start the flow^~^ . 
Summary : 

It is recommended to use the downstream diaphragm location, 
mainly because the starting time in the subsonic portion of the nozzle 
is very short ( - 1 msec) and the troubles of having diaphragm particles 
in the flow are eliminated. For supersonic experiments, it might be an 
advantage to have the diaphragm located upstream. The diaphragm holder 
unit can be designed so that there will be a possibility to mount it 
either downstream or upstream. 

4. Charge tube pressures and temperatures 

The pressures and temperatures, and , in the charge tube, prior 
to diaphragm rupture, as a function of pressure and temperature T^, 
respectively, where the Mach number of the flow is , is given by^^ : 
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(1) 



Y - 1 2 

!«!i ,!i V !2 1 '^~ m 3 

T 4 f st T 4 T 3 T st'(l + i^M,) 2 






. a . <_v - 1 

P 4 T 4 



( 2 ) 



2 1 

For Y = 1.4, M = 0.2 + — = 0.93 _ , J 

J 4 1 + Y - M 2 

T* 

1 



For y = 1.67, M = 0.2 + — = 0.95 
J T 4 






The corresponding values of P^/^ can ca l cu l ate< ^ from Eq. (2) . 

P, as a function of P n for different values of the M_ is shown 
4 2 2 

in Fig. 5 and Fig. 6 for y = 1.4 and y = 1.67 respectively. The 



values of T for y = 1.4 and y = 1.67 and different values of are 



given in Table V assuming T^ = 300 k: 



M 2 


0 


0.5 


1 


1.5 


2 


Y = 1.4 


280 


266 


233 


193 


155 


Y = 1.67 


285 


263 


213 


162 


122 



Table V 



The pressure, P^, in the part of the system downstream the diaphragm 
should be of the order of P^/100 (to get best start times without a 
diffuser) ; i.e. P^ ~ 1 torr (see Sec. 5) . 

5. Diffuaar, dump tank and valve 

The starting times, as discussed before, become longer with longer 
diffusers, booh for downstream and upstream diaphragm locations. 
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For a fixed charge tube pressure , and back pressure being in- 
creased gradually, the flow at the diffuser end ceases to start for a 
certain value of while further increasing P^, the supersonic flow at 
the nozzle exit barely starts. Beyond these critical conditions, the 
flow at the nozzle exit becomes unstable. The limiting experimental 

(15) 

downstream pressure radios, P^/P^, were measured by Merkli and Abuaf 
both for upstream and downstream diaphragm locations, as a function of 
diffuser length L/D. 

It was found that (P^/P^) cr for diffuser starting is -0.11 inde- 
pendent on L/D. 

(P^/P^) for nozzle supersonic flow increases with L/D and reaches 

a maximum at about L/D - 6 where (P,/? J =0.22. 

1 4 cr 

That means that within the available time, the pressure in the dump 
tank should not increase above 0.1 x P^; this limiting pressure defines 
the volume of the dump tank. 

The rate of pressure rise inside the tank: The continuity dictates 

* p d 

that in the dump tank V , t + m _ =0 where V, is the tank's volume, p. 

^ d d t d d d 

is the density and m^ is the total mass flow into the tank. By differen 
tiating the equation of state, ?^ = R T^, with respect to time and 
using the continuity equation, we get 



0(r - fji . 

3t RT d 3t “ V d m d 



A d = *3. = P 3 A 3 V 3 = 



P 3 A 3 M 3 a 3 
RT„ 



T^S^A^M^a- 



A _ uT3 3 3 

AP , = — — T 

d TV 

3d 
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Where x is the available testing time. Since initially P d = P^ << 

0.1P„ we may write AP, = P.; from the last equation: 

3 da 

T p 

V' . = (— ) (— : )A M <. 2 x 10 x n x 8 2 x0.2 x 33000 x 20 x 10~ 3 
d T 3 P d 3 

V - 5.3 x 10 5 an 3 
a 

It is suggested to install a valv.e upstream the diapragm, as 
shown in Fig. 3. After each run, before the ruptured diaphragm is 
taken off, this valve will be closed and the gas in the system will 
be saved. 

A more efficient way of saving gas will be to use a compressor 
to compress the gas in the dump tank and recirculate it back to the 
charge tube. The compressor should be vacuum tight. A simple and 
inexpensive compressor is suggested in Fig. 7. 
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IV. INSTRUMENTATION 



Pressure gauges 

0-5 atm. range for measuring . 

0-100 torr ^ 

0-760 torr Pressure gauges for gas mixtures preparation . 

0-2 atm j 

0-100 torr for measuring P^ 

For measuring P^ and P^ during the experiment itself, after 
diaphragm is broken, we need piezoelectric pressure transducers with 
time response of the order of tens microseconds, in the range of 
0-2 atm. 



Velocity profiles and turbulence spectrum 

For transient flow speed and turbulence measurements, a laser- 
doppler anemometer or a hot-wire anemometer can be used. The avail- 
able steady flow time, 10-20 msec, might be too short for appropriate 
averaging the turbulence spectrum, mainly in the low frequency regime. 
Therefore it will be necessary to run the flow in a continuous mode (with 
compressed air) for turbulence spectrum measurements. The transient flow 
fluctuations may not be averageable but their ability to stabilize the 
discharge may surpass that of steady flow in which case other measurement 
forms should be sought. 

Voltage current measurements ; probes 

The voltage and current have to be recorded during the experiment. 
The simplest way of doing it would be by using an oscilloscope. 

The probes used for measuring plasma parameters must be provided 
with automatic bias sweep, so that the characteristics can be 
accomplished within the flow time. 
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Lasing diagnostics 



A probe laser for s.s.g. measurements. 

A fast response detector (microseconds range) . 
Power meter. 
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VACUUM PUMP 
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Fig 2. TRANSIENT FLOW SYSTEM SCHEMATIC. 



CROSS SECTION A-A 
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Fig 3. SUGGESTED TURBULENCE GENERATING PLATES MOUNTED CLOSE 
TO ELECTRODE TIP. 



DIAPHRAGM 





FLOW 




Fig. 4(a) APPROXIMATE WAVE DIAGRAM FOR THE FLOW IN A 
LUDWIEG TUBE. UPSTREAM DIAPHRAGM LOCA- 

TION. 
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DIAPHRAGM 




Fig. 4(b) APPROXIMATE WAVE DIAGRAM FOR THE FLOW IN A 

LUDWIEG TUBE.- DOWNSTREAM DIAPHRAGM LOCATION. 
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VACUUM TIGHT COMPRESSOR. 
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